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Recent work in two brane Randall-Sundrum models has shown that the radion can drive paramet¬ 
ric amplification of particles on tfie TeV brane leading to a stage of cosmology similar to preheating. 
We consider radion induced preheating as a mechanism for electroweak scale baryogenesis and es¬ 
timate the baryon to entropy ratio for some regions of parameter space. We also predict an upper 
bound for the radion mass, which makes this baryogenesis mechanism testable. Finally, we verify 
some assumptions with numerical calculations of the H-l dimensional Abelian Higgs model. 


I. INTRODUCTION 

The high degree of symmetry between matter and anti¬ 
matter in our theories has long seemed in contradic¬ 
tion with the overwhelming presence of matter (baryons) 
around us. Electroweak baryogenesis 0 provides a 
testable explanation of the baryons within the context 
of the standard models of particle physics and cosmology 
but suffers from the two difficulties of providing enough 
CP violation and providing enough departure from equi¬ 
librium. In the usual picture transitions between vacua 
of the electroweak theory correspond, via an anomaly, to 
changes in baryon plus lepton number. If the electroweak 
phase transition in the early universe is strongly first or¬ 
der then the presence of CP violation will cause sphaleron 
transitions near bubble walls to favor formation of mat¬ 
ter over anti-matter. While the minimal supersymmetric 
standard model may provide the first order phase transi¬ 
tion and CP violation needed for this picture of baryoge¬ 
nesis, a radically different picture is possible within the 
context of the Randall-Sundrum model. 

As with supersymmetry, the two-brane Randall- 
Sundrum model 0 (R-S) provides a framework for ad¬ 
dressing the hierarchy problem. The model consists of a 
slice of five dimensional anti-de Sitter space (AdS^) trun¬ 
cated by two 3-branes with the standard model particles 
trapped on one of the branes. The gravitational warping 
between the branes generates hierarchies of mass scales. 
With the stabilization of the brane separation provided 
by Goldberger and Wise 0, mass scales are naturally 
stabilized providing a solution to the hierarchy prob¬ 
lem along with an acceptable low-energy phenomenology. 
The separation between the branes corresponds to a par¬ 
ticle called the radion which is expected to be the lightest 
new degree of freedom. 

Initially, it might seem difficult to explain the baryon 
asymmetry within the context of RS since RS physics 
becomes important just above the electroweak scale. For 
example, it is not possible to use existing extensions of 
standard model physics to explain baryogenesis within 
an effective theory that integrates out RS physics. They 


are at the same scale. We must, therefore, re-examine 
the baryogenesis question within the context of the RS 
model. Fortunately, the radion provides a new route 
for departure from equilibrium; previous work has shown 
that homogeneous oscillations of the radion lead to expo¬ 
nential growth of scalar particles on the standard model 
brane Q. The produced particles do not have a ther¬ 
mal distribution, but lie in a small region of momentum 
space. These are features common to parametric ampli¬ 
fication, or preheating [^ , as it occurs in the inflationary 
context. 

We use parametric amplification of the standard model 
Higgs field to drive the field over the potential barrier, 
allowing for the creation of winding modes and, conse¬ 
quently vacuum transitions. This parametric amplifica¬ 
tion, rather than bubbles from a first order phase transi¬ 
tion, generates the departure from equilibrium. With an 
additional source of CP violation, the decay of the Higgs 
windings will favor the production of baryons over anti- 
baryons. However, there are constraints on the mass of 
the radion coming from the size of windings which give 
rise to baryon number violating transitions. Generically, 
if the radion is too massive, then the parametric amplifi¬ 
cation gives windings which are too small for baryogen¬ 
esis. The radion should not be much more massive than 
three times the mass of the Higgs boson. Furthermore, 
because the radion is a gravitational mode in the full the¬ 
ory, the form of the radion-Higgs coupling is determined 
with only one free parameter (we consider only minimal 
coupling of the Higgs). 

The next sections review relevant aspects of the stan¬ 
dard model baryon violation and the RS model. We then 
discuss radion induced preheating and consequences for 
baryogenesis. We explain some constraints on the model 
and estimate the baryon asymmetry for a simple case. 
The final section shows numerical results which verify 
some assumptions and help develop our intuition. 


II. BARYOGENESIS IN THE STANDARD 
MODEL 


- As alluded to above, non-perturbative effects in the 

electroweak sector of the standard model violate baryon 
‘Electronic address: mmartin®cmu.edu number. An anomaly links the change in baryon and 
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lepton number to the change in Chern-Simons number: 

AB = AL = nfANcs, ( 1 ) 

where n/ is the number of families of matter fermions. 
In the absence of fermions the vacua of this theory are 
degenerate and they are labeled by the winding num¬ 
ber of the Higgs field. In any vacuum state the Chern- 
Simons number is equal to the Higgs winding number. 
The transition from one vacuum state to another must 
therefore involve a change in both the Higgs and gauge 
fields so that both the Higgs winding and the Chern- 
Simons number remain equal. See the reviews in P] for 
more details. This baryon violating transition between 
vacua passes through a field configuration known as the 
sphaleron, which is a saddle point of the effective po¬ 
tential. Henceforth, we refer to the sphaleron mediated 
transitions between vacua as the sphaleron transition. 

We can now understand that an initial configuration 
with the Higgs winding differing from the Chern-Simons 
number by one unit must relax to a vacuum in one of 
two ways: either the Higgs winding will change, which 
does not violate baryon number; or the gauge field evo¬ 
lution will change the Chern-Simons number, which does 
violate baryon number. The method by which the re¬ 
laxation occurs is dependent upon the length scale of the 
winding. In general a Higgs winding will collapse and the 
time scale for the collapse is simply the length scale [^. 
However, the gauge fields will evolve to cancel gradient 
energy in a Higgs winding with time scale loosely set by 
the gauge boson mass. So a sufficiently small winding 
will collapse and the Higgs winding number will change, 
while a sufficiently large winding will be canceled by the 
gauge field thus changing baryon number. 

Some critical length scale separates these two regimes. 
If there is initially an equal density of windings and anti¬ 
windings on large scales and no CP violation, then there 
will be equal production of baryons and anti-baryons. 
According to the picture of Turok and Zadrozny [^, 
CP violation will affect the critical length differently for 
windings and anti-windings. An anti-winding near the 
critical scale will preferentially unwind through Higgs 
field evolution, while a winding at the same scale will 
unwind by gauge field evolution, thus favoring produc¬ 
tion of baryons. The critical length is about ~ 

(30 GeV)~^ if the gauge boson mass is given by the zero 
temperature vacuum. The realistic picture is more com¬ 
plicated 0 with the length scale being a function of 
other parameters which characterize a winding, but there 
should still be a range of critical length scales separating 
Higgs field unwinding from baryon number changing. 

In the high temperature, symmetry restored vacuum, 
simple dimensional analysis gives the sphaleron length 
scale as being the magnetic screening length, (awT)~^, 
which could be as large as about (3GeV)~^. For the 
remainder of the paper we work with the length scale 
Lph = (20GeV)~^. We use this length scale to constrain 
the radion mass since a massive radion will give higher 


frequency excitations to the Higgs field leading to wind¬ 
ings which are too small for baryogenesis. 

At temperatures much below the electroweak scale, 
sphaleron transitions are exponentially suppressed by the 
energy scale of the potential barrier between electroweak 
vacua, Esph ~ lOTeV: 

( 2 ) 

For temperatures near the electroweak scale and as low as 
150 GeH Hi, the electroweak symmetry is restored and 
the sphaleron transitions are unsuppressed.^ Thermal 
fluctuations easily form windings of the Higgs or gauge 
fields and on dimensional grounds 

Tsph ~ iawT)\ (3) 

It might seem then that a temperature above 150 GeV, 
along with CP violation, is sufficient for the production of 
baryons. However, the thermal average of baryon number 
is zero 0 , meaning that baryogenesis must happen dur¬ 
ing a departure from thermal equilibrium. Furthermore, 
if equilibrium is restored while sphaleron transitions are 
unsuppressed, then any existing baryon asymmetry will 
be destroyed. This process is called thermal washout, 
and it will constrain the total energy in our model as 
explained in section nvi 

III. THE RANDALL-SUNDRUM SCENARIO 

The background geometry for the two brane Randall- 
Sundrum model (RS) is given by the metric 

+ dz^, (4) 

with branes truncating the space at z = zq and z = zi. 
This geometry is a solution to Einstein’s equations for 
negative bulk cosmological constant and nonzero brane 
tensions which scale with k. For naturalness k is taken to 
be near the fundamental Planck scale, M 5 , but somewhat 
smaller to ensure gravitational perturbativity. The warp¬ 
ing of the space red shifts mass scales on the brane at zi 
relative to the scale at zq, allowing for the electroweak 
theory trapped on the brane at zi to have an energy scale 
much less than scales on the brane at zq. The brane at zi 
is therefore often called the TeV brane. The separation 
between the branes corresponds to a modulus known as 
the radion. By fixing the radion in a natural way, the 
electroweak hierarchy problem is solved. 

The form for the radion which does not mix with the 
four dimensional graviton 0 is 


^ Some authors Q| cite 100 GeV as the temperature above which 
sphaler on transitions are unsuppressed. We choose to follow ref¬ 
erence since a larger temperature is more conservative with 
respect to placing bounds on the radion mass. 
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+ {l + 2kf{x)e^'^^)^dz^, (5) 

where /(x) is the radion. Goldberger and Wise 0 found 
that a scalar field in the bulk which has nonzero vacuum 
expectation values on the branes can naturally determine 
the brane separation. A competition between the kinetic 
and potential terms in the five dimensional action for the 
scalar gives a potential to the radion in the four dimen¬ 
sional effective theory. In other words, after integrating 
out the bulk degrees of freedom, an effective theory for 
the radion remains 

= -^df,Fix)df^Fix) - ^mlF{xf + 0{F^). ( 6 ) 

We have rescaled fix) to have the properly normalized 
kinetic term, 

Fix) = ^Jl2kMle^^^'^-^°'>fix) (7) 

and the value for the radion mass, niF, depends on the 
details of the bulk scalar. The scale of physics o n the TeV 
brane is set by the parameter A = 3M| /ke~^^^'^-^o) _ 
We will only consider small oscillations of the radion. 
Fix) ^ A, to avoid sensitivity to the higher order terms 
in the radion effective action. 

Now consider a scalar field on the TeV brane. We will 
see how the warping red shifts mass scales and derive the 
coupling of the radion to the scalar needed for the next 
section. We set zq = 0 for simplicity. The action on the 
brane is 

S' = - i y d'^x^/^ + rmfcj?] , (8) 


with the induced metric on the brane given by 

ds^ = g^^ix)dx^dx'^ = ^'^ij^^dx^dx''. 

( 9 ) 

By rescaling (/)(x) —> e^^^cfix) to give the kinetic term 
canonical normalization, the action takes the form: 


S = -- 


1 


d*x 






■ ( 10 ) 


Notice that the scalar field mass is suppressed by a “warp 
factor”, e~^^^. We define the effective scalar mass = 
e-kzim. In terms of the correctly normalized radion, the 
action is: 


IV. PREHEATING AND BARYOGENESIS 


As we discussed in previous work 0, there may have 
been a time in the early universe when most of the en¬ 
ergy density was tied up in coherent homogeneous oscilla¬ 
tions of the radion. These oscillations might be expected 
for several reasons. First, an inflation-like mechanism is 
needed to solve the horizon problem within the RS con¬ 
text, and this mechanism almost certainly will involve 
physics of the extra dimension. In a minimal model the 
only way to reheat the TeV brane is through the sta¬ 
bilization mechanism, which would involve oscillations of 
the stabilizing field and of the radion. We take the poten¬ 
tial of the stabilizing field to be stiff and study oscillations 
of the radion. 

Alternatively, at high enough temperatures the space- 
time is unstable to formation of a bulk horizon, which 
cuts off the space and removes the TeV brane. As the 
universe cools, a phase transition occurs |3, generating 
the TeV brane generically shifted from its equilibrium 
position. The radion will oscillate coherently as the brane 
settles in to the equilibrium position. In this scenario we 
would not expect a cold TeV brane, although more de¬ 
tails of this phase transition still need to be understood. 

In a third picture, the stabilization mechanism may 
provide a false vacuum for the radion with the TeV brane 
at infinity ^ 3 . Early universe dynamics could set the 
brane in this false vacuum with a first order phase transi¬ 
tion bringing the brane to the global minimum. Bubbles 
of the brane in the true vacuum would not be expected to 
lie exactly at the bottom of the potential, but oscillate 
around the minimum. In any of these three scenarios, 
higher energy physics leads to radion oscillations. 

We take the radion to be a background of the form 

Fit) = FoAcosimpt), ( 12 ) 


and decompose the scalar field in Fourier modes: 
d^k 


(j)it,x) = 


(2^)3/2 I 


^<fkit)' 


-ik-: 


Ulit)e^'^"- 


(13) 

We then calculate the equation of motion for the scalar 
modes from the action Hill : 

-I- Fomp sinimpt)^^ + Jlit)(j)k = 0, (14) 

dV dt 


S = -- 




where w^(t) = 

|fc|2 -hTO2g-Focos(miri) 

/ 


vacuum state for all modes, 


We start with a 




( 11 ) 


This is the action used to study parametric amplification 
of scalars on the TeV brane |^. Notice that the univer¬ 
sality of gravitational couplings limits the radion-scalar 
coupling to this form. In the next section we review this 
parametric amplification and show how baryogenesis oc¬ 
curs. 


ddh 

^(0) = -*^fc(0)</>,(0), <^,(0) = ^==, (15) 

dt y/2uJkiO) 

and find exponential growth of some modes. For the 
specific case Fq = 0.5, the modes unstable to growth are 
shown in the dark region of figure ^ This is parametric 
resonance, or preheating. 
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FIG. 1: The instability band for the massive scalar when 
Fo = 0.5. The dark region shows where there is exponential 
growth of scalar particles. Both mp and k = \k\ are expressed 
in units of the scalar mass. 


It should be emphasized that the spectrum of parti¬ 
cles produced through parametric amplification is non- 
thermal and the time scale for growth of particles is given 
by the radion mass. To show this we calculate the expec¬ 
tation value of the number density operator, {Nk{t)) (see 
the appendix for details). In figure |21 we plot \n{Nk{t)). 
The modes in momentum space which grow are those 
expected by comparing with figure ^ We expect that 
when backreaction of the scalar on the radion is taken 
into account, the scalar will drain energy off the radion 
and damp the amplitude of radion oscillations. 

We now want to think about placing the Standard 
Model on the TeV brane. Again assuming homogeneity 
and an initial state near the zero temperature vacuum, 
we would expect Higgs particle production as explained 
above, where is now the Higgs mass. Of course, the 
action for the Higgs field should involve the self-coupling 
and the coupling to gauge fields. However, leaving these 
terms out should be a good approximation at the begin¬ 
ning when only the minimum of the potential is being 
sampled by the Higgs vacuum expectation value. The 
choice of initial conditions is generalized in the numeri¬ 
cal work. We discuss production of fermions and gauge 
fields below. 

The essence of the idea is that the preheating of the 
Higgs field will produce configurations with winding, in 
a manner analogous to the production of topological de¬ 
fects following inflationary preheating. These windings 
will then decay. Some decays will violate baryon number, 
some will not, and the CP violation will favor production 
of a small number of net baryons. A similar idea of in- 



FIG. 2: Log of the number density operator. In Nk{t), for 
niF = dm^i, and Fo = 0.5. Note that for momenta within the 
instability band of hgure0 the number of particles produced 
grows exponentially. All units are appropriate powers of the 
scalar mass, nir/,. 

flationary preheatin g le ading to baryogenesis has been 
proposed by others [Ifl 1111 ■ However, radion preheat¬ 
ing occurs naturally at the proper energy scale, the elec- 
troweak scale, and has little freedom in parameter space 
due to the universality of radion couplings. These con¬ 
cepts are significant: there is no fine tuning needed, as in 
the inflationary case, to make radion preheating happen 
at the correct energy scale; and our model will be tested 
and may be excluded at the next generation of colliders. 

A. Constraints 

Before estimating the amount of baryogenesis, we 
would like to discuss some constraints on the radion. 
First, the final reheat temperature should be sufhciently 
below the electroweak phase transition to avoid thermal 
washout of the baryon asymmetry. If, once thermal equi¬ 
librium is reached, we require 

^reheat ^ 150 GeP, (16) 

then we must constrain the total energy density: 

2 

p = < (350 GeVf. (17) 

Because the expansion rate of the universe, H ~ ^ 

10“^®GeH, is much less than other relevant scales in the 
problem, we may ignore the expansion. Then energy con¬ 
servation restricts the initial radion configuration: 

PF = <{?,mGeVf. (18) 

Another constraint comes from the location of the in¬ 
stability band. As explained in section^ the length scale 
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FIG. 3: The instability band for the massive scalar when 
k = 0.16 m^. This corresponds to A: = for a Higgs mass 
of 125 GeV. The dark region shows where there is exponential 
growth of particles. The radion mass is expressed in units of 
the Higgs mass. Above the line the system contains enough 
energy to washout baryons after thermalization. 

for sphaleron configurations, Isph, is about (20Gey)“^. 
The instability band therefore needs to encompass this 
value of momentum in order to produce sphaleron config¬ 
urations. For the case considered in figure^ this criteria 
is only satisfied for a radion mass of 1.5 < < 3. 

More generally, we can understand this constraint by 
plotting the instability band as a function of Fq and mp 
for fixed k = as in figure 13 

The region above the line is excluded by the en¬ 
ergy density constraint. This assumes a Higgs mass of 
125 GeM, to set a scale, and a fairly low value for the 
cutoff, A = 750Gel3. These values are used in the re¬ 
mainder of the paper. Equation m shows that this line 
moves down (excluding more of the preheating band) as 
the lower limit of the cutoff increases or if the Higgs mass 
becomes heavier. Already we can see that if the radion 
mass is much above three times the Higgs mass, then 
there can be no preheating without thermal washout of 
baryons. For a light radion we need large values of Fq 
for preheating. To trust a calculation in this regime re¬ 
quires knowledge of the complete radion potential. This 
involves details of the stabilizing mechanism. 

While 0{F^) and higher terms in the radion action 
will alter the form of radion oscillations and make a pre¬ 
cise large Fq calculation of the baryon to entropy ratio 
difficult, in practice, these terms do not prohibit para¬ 
metric amplification of the Higgs field. A simple numer¬ 
ical calculation similar to that shown in figure 13 demon¬ 
strates that arbitrary 0{A) and 0(1) coefficients for F^ 


and F^ terms respectively preserve exponential growth 
of the Higgs fluctuations. The significant concern for the 
purposes of this paper is that higher terms would modify 
the location of the instability band without a comparable 
modification of the energy density constraint, thus alter¬ 
ing the upper limit on the radion mass. While we can’t 
rule out this possibility, it does not seem to be a generic 
feature of these operators. 

B. Evolution 

We may estimate when sufficient preheating has oc¬ 
curred to form Higgs windings by calculating the two 
point correlation function of the Higgs fluctuations. We 
work in a toy model with no gauge fields. Shifting the 
Higgs field, 

</' = F + X, (19) 

the equation of motion for the vacuum expectation value 
(VEV) of the Higgs in the Hartree approximation m 
becomes: 

n(p + X{ip'^-v'^+ 3{x^))(p = 0. (20) 

In the electroweak theory windings will be able to form 
when the Higgs VEV crosses the potential barrier. In 
this toy theory, ip crosses the barrier when this two point 
function is of the order of the zero temperature vacuum 
expectation value: 

(x^)~uV3. (21) 

It is not clear that the system will necessarily behave as if 
the symmetry of the vacuum for p has been restored, but 
we do not need to concern ourselves with this issue. We 
only need the fluctuations, (x^), to kick the field across 
the potential barrier for winding formation. 

In figure 0] we plot (x^)(t) — (x^)(0) where 

(x^)(t) = ^^^lxfc(i)r (22) 

The equation of motion for the Eourier modes, Xk, is 
given by equation la where (jik is replaced by Xk and 
m,j, is now the Higgs mass. We start again with vacuum 
initial conditions as in equation Erom the plot it 

is clear that the winding formation condition, equation 
is satisfied after a few radion oscillations. Of course, 
at this point the approximation of linearizing the equa¬ 
tions has broken down. We are simply making the point 
here that windings should form on a time scale set by 
the radion mass. Because (y^) is cutoff dependent even 
in the initial vacuum state, we plot the change in this 
quantity. 

There are several possible times scales relevant to this 
problem. Eirst, the Hubble scale, is often important 
in electroweak baryogenesis since the Hubble expansion 
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FIG. 4: The two point function for the Higgs fluctuations 
with mp = 2.25 and Fq = 0.5 with units in powers of the 
scalar mass. The horizontal line is n^/3 when ra^ = 125 GeV. 

cools the universe and causes the sphaleron transitions 
to freeze out. Here that role is played by the thermaliza- 
tion time scale and the Hubble parameter is irrelevant. 
The sphaleron transitions freeze out as the energy in the 
Higgs field is redistributed to all the other degrees of 
freedom in the standard model. The total decay rate of 
the Higgs boson is near 10“^ GeV 0 so we assume this 
determines the thermalization time scale. We take the 
sphaleron time scale to be the length scale, Isph, which is 
much shorter than the thermalization scale. Finally, the 
fourth time scale is that of preheating. As we have seen 
above, the preheating scale is given by the radion mass. 

This picture is not significantly altered by the pres¬ 
ence of matter fermions or gauge fields. Pauli blocking 
prevents the exponential growth of fermions, and the de¬ 
cay time scale of the Higgs boson is much longer than 
the preheating scale. Fermions can safely be ignored. 
The radion also couples to the gauge fields, but not to 
the kinetic term. Without the derivative coupling that 
was present between the radion and the Higgs in equa¬ 
tion m, there is much less parametric amplification. 
The full sphaleron transition naturally involves gauge 
fields, but here we focus on generation of the Higgs wind¬ 
ings only and leave more detailed study for future work. 


C. An Estimate 

For generic parameters inside the preheating band of 
figure 13 the radion will continue to produce windings 
of the appropriate length scale for much longer than the 
sphaleron time scale. This winding production time will 
depend on the total radion energy and on the thermaliza¬ 
tion time scale for the system. For example, in hgure^l 
only about 1/10 of the energy of the radion has gone in 
to the Higgs field by the time windings form. We ex¬ 
pect baryon production to be quite efficient here, since 
sphaleron transitions may take place throughout the en¬ 
tire brane during preheating. To calculate the baryon 
asymmetry, though, will require a deeper understanding 
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FIG. 5: The mp = 2.5 plane of the instability band. The 
momentum, k, is in units of Higgs mass. 


of the sphaleron transition rate in this far from equilib¬ 
rium state. It will also require the full SU (2) Higgs and 
gauge fields, a better understanding of the thermalization 
time scale, and, of course, a source of CP violation. There 
has been work on this calculation within the context of 
baryogenesis following inflation. Some efforts have ap¬ 
plied ideas from thermal sphaleron transitionsHlJ , while 
others have applied lattice calculations [H ll 5|| . These 
steps are important, but difficulties may remain |Hi|. 

Nevertheless, for some specific choices of Fq and mp, 
it may be possible to estimate the winding density of the 
Higgs field and acquire insight on the baryon asymmetry 
in a manner which is independent of a new CP violating 
mechanism. If the radion preheating generates windings 
and shuts off on a time scale shorter than the sphaleron 
time, then the Higgs winding density can be estimated 
via the Kibble mechanism. To see how this might hap¬ 
pen, we plot the instability band as a function of Fq and 
momentum in figure [3 for radion mass mp = 2.5 m^. 
Considering the backreaction of the Higgs field on the 
radion, we expect a smaller Fq to be relevant later as the 
radion loses energy to the Higgs field. 

So, if the system starts near Fq = 0.3, then sphaleron 
sized windings may occur for a time Isph- As Fq de¬ 
creases, much of the remaining energy of the radion will 
be dumped to higher momentum modes centered around 
k = 0.75 m^. Then, much as happens in defect formation 
through inflationary preheating, we expect that there will 
be order one winding per region of space. In this pro¬ 
cess the quantum fluctuations are being amplified, and 
unless there is reason to expect long range correlations, 
then causally disconnected regions should wind in differ¬ 
ent directions. 
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We may now make our estimate. We calculate the 
entropy present: 


27r2 

s = ~^9*sT , (23) 

where the temperature after thermalization is given by 
the energy density for Fq = 0.3, and mp = along 

with equation 113: T Ri 0.75 m^. We parameterize the 
CP violation such that the baryon density is the product 
of the winding density and ecp- 


B = UwCc 


Then the baryon to entropy ratio becomes 




cp 


S 


lO-'^ecp 


(24) 


(25) 


The measured value is ry ^ 10“^°. Note that the strong 
dependence of the winding density to the sphaleron time 
scale makes this estimate imprecise. 

It is common in four dimensional scenarios to postulate 
that a nonrenormalizable CP violating term arises from 
a higher energy theory. The same may be possible here, 
but any higher energy theories need to live in the bulk. It 
would be interesting to investigate the details of bulk CP 
violation and the consequences for baryogenesis on the 
brane. It may also be possible to add the CP violation 
on the brane at electroweak energies, as in a two Higgs 
doublet model. 


V. A NUMERICAL MODEL 

To see some details of this process borne out more ex¬ 
plicitly, we now turn to a numerical implementation of 
the 1-1-1 dimensional Abelian Higgs model. This model 
has features in common with the baryon violation in the 
standard model, as we will see below, and has therefore 
often been used in the past to study baryogenesis 0. 
The Abelian Higgs model allows us to go beyond the 
linearized equations and to include backreaction on the 
radion without the technical difficulties of solving a non- 
Abelian model numerically. We can study the forma¬ 
tion of winding modes and confirm our intuition that the 
length scale of the windings is given by the location of 
the instability band. 

The action for our toy theory is 

S = J d^x(^- 

(26) 

In the vacuum the scalar sits at the bottom of the poten¬ 
tial well 


(p = ve 


i9(t,x) 


(27) 


and is pure gauge 

Ap = ^dp0{t,x). (28) 

This theory has a winding number which labels the vacua 

NH{t) = ^Jdxd,9, (29) 

and transitions between vacua require (p to pass through 
zero. In addition, the Chern-Simons number is propor¬ 
tional to the Higgs winding number in the vacuum state: 

Ncs = j dxA, = fNH. (30) 

Thus the vacuum structure of this theory resembles that 
of the electroweak sector of the standard model. With the 
addition of fermions there is a conserved chiral current, 
which is anomalously broken in a way that is related to 
the change in Chern-Simons number [2, A{Nfi — Np) = 
ANcs- 

Because the radion, which stabilizes a slice of AdSs, 
does not couple as the 5 dimensional radion, we intro¬ 
duce a “radion” designed to mimic the 3-1-1 dimensional 
preheating. As mentioned above, we are primarily in¬ 
terested in the early time behavior and specihcally the 
formation of the Higgs windings modes. Therefore, we 
restrict our attention to the radion-Higgs sector: 

S = J d^x(^—\e~^/^'q^^''dp(p'^du(p 

_Ag-2F/A _ ^2)2^^ (31) 

We start the Higgs field in a vacuum state as be¬ 
fore. We then briefly evolve the linearized quantum equa¬ 
tions da- After particle production has begun, but be¬ 
fore equation m is satisfied, we Fourier transform to 
configuration space and evolve the full non-linear clas¬ 
sical equations of motion. In other words, the quantum 
field with a few particles is serving as the initial conhgura- 
tion of the classical held for our classical evolution.^ We 
checked that results are not sensitive to the time when 
the transition from quantum to classical evolution oc¬ 
curs. Because of the exponential growth in the preheat¬ 
ing band, results are also relatively insensitive to random 
changes in the initial conditions, equation 113, by fac¬ 
tors of 4 or more. With such large initial fluctuations, 
there are more windings at early times, but after several 
oscillations of the radion there is not a signihcant dif¬ 
ference. However, an initially hot brane, or a brane with 
any non-trivial initial energy density, does make a change 
to the reheat constraint, equation m since it is the to¬ 
tal energy density which is constrained by equation 03. 
So initial energy density on the brane will move the line 


^ After our work was completed, reference Hi appeared using the 
same method of setting initial conditions for the classical field. 
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FIG. 6: Phase of Higgs field shortly after windings begin to 
form. The length, x, is in units of inverse Higgs mass. 

in figure 13 down, thus excluding more of the preheating 
band. 

In figure El we show the phase of the Higgs field af¬ 
ter windings have begun to form. We want to know the 
length scales of the windings. At early times, as in fig¬ 
ure El windings are easy to identify since the fluctua¬ 
tions in phase have not grown too large. We define the 
length of winding as being the minimum distance be¬ 
tween two points separated in phase by 27r. In figure [7| 
we show the number of windings as a function of length 
scale for two cases with Fq = 0.5, as in figure Q] The 
curve peaked at smaller scales corresponds to the more 
massive radion, which has a preheating band at larger 
values of momentum. This is what we expect. For the 
case with longer length windings, mr = 2.0 m^, there 
are about 1200 windings. This corresponds to about one 
winding in every 10 causal volumes. The other scenario, 
niF = i-Om^, preheats more efficiently, in part because 
the radion oscillates more quickly leading to about one 
winding in every 2 causal volumes. In both cases the 
winding density is increasing quickly with time. 

It would be interesting to use numerical methods on 
the classical 3-1-1 dimensional SU{2) theory to track the 
Chern-Simons and Higgs winding numbers. The depen¬ 
dence of these quantities on the radion mass should un¬ 
ambiguously confirm or refute the model. In addition, 
this calculation would provide a more concrete estimate 
of the range of baryon to entropy ratios possible. 

VI. CONCLUSION 

We have shown that baryogenesis can occur in the 
two-brane Randall-Sundrum model with radion preheat¬ 
ing providing the departure from equilibrium. Another 
source of CP violation is needed, and as mentioned, inves¬ 
tigations of bulk CP violation may prove interesting. One 
of the successes of this model is that preheating occurs 
naturally at the energy scale for sphaleron transitions to 
generate the asymmetry, but not wash it out. 

Furthermore, the existence of a baryon asymmetry in 


FIG. 7: Number and length scale of windings a few time units 
after winding formation in units of Both use Fb = 0.5, 

but have different radion mass: mp = 4.0 on the left 
and mp = 2.0 on the right. The length scale relevant for 
baryogenesis, laph, is in the range of The total length 

for the simulation was 10® 

our universe implies an upper bound on the radion mass. 
For a sufficiently large radion mass, the preheating band 
does not extend down to sphaleron length scales, unless 
the radion energy exceeds the thermal washout limit. We 
conclude that the radion mass should not be much larger 
than three times the Higgs mass, if the baryon asymme¬ 
try is to arise from this model. 

APPENDIX: NUMBER DENSITY OPERATOR 

We calculate here the expectation value of the number 
density operator, (0|aj^(t)afc(t)|0), for a massive scalar 
field coupled to the radion, i.e. for action CD- First, let 
us change the notation slightly to ease this calculation. 
We decompose the scalar: 

/ d^k 

(A.l) 

where 

^k{t) = Uk{t)dk + (A.2) 

We decompose the conjugate momentum, Tr{t,x), in the 
analogous manner and find 

^kit) = {Uk{t)ak + Ulkit)d^-k) ■ (^-3) 

We have assumed that the radion is independent of the 
spatial coordinate in anticipation of homogeneous oscil¬ 
lations. 

This decomposition is in terms of operators which an¬ 
nihilate the initial vacuum, afc|0) = 0. We seek an ex¬ 
pression for time dependent operators, dk{t), in terms of 
the initial time operators so that we may calculate the 
number of particles relative to the initial vacuum. The 
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time dependent (Heisenberg picture) operators must sat¬ 
isfy the equation of motion 

= [ak{t),H{t)], (A.4) 


Finally, using the standard equal time commutation rela¬ 
tions for the field, •7r^(t)] = i{2Tr)^S^{k+^ we may 

evaluate the equation of motion for dk{t). The result is 
simple: 


where H {t) is the Hamiltonian. The Hamiltonian can be 
written in terms of (j)k{t) and TTk{t)'. 


lg{t) = 0=lh{ty (A.9) 


H(i) = 




2 J (27r)3/2 L 


e^(*)/A7rfc(t)7r_fe(t) 


+e 




^ki't)(l^k{t)(l)-k(t) . (A.5) 


We now use the form for the oscillating radion given in 
the text in equation 03 and evaluate the number oper¬ 
ator from the ansatz 1A.6II : 


As in the main text, = k'^ + m^e We may 

look for a solution by making the ansatz: 

afc(i) = g{t)%k{t) + h(t)(j3k{t). (A.6) 

This guess is motivated by the form of a Bogolubov trans¬ 
formation and the knowledge that 7 !^ {t) and (j)k (t) contain 
all of the degrees of freedom of the system. 

From the initial conditions, dk{t = 0) = dk, and 

^(0) = -7..(0)f/.(0), UkiO) = -^^, (A.7) 

as in equation 03 , we find the initial values: 

5(0) = lUm h{0) = -7e-^(°)/^f7,*(0). (A.8) 


(Nkit)) = (0|a[(t)afc(t)|0) = 


o — Fo 


2uJk(0) L 


32 f’o[l—cos(mFt)] 


\Uk{t)\^+ioimuk{t)\ 


le 


— Fq cos(mFt) p 


u^mkit)-u:it)Ukit) .(A.io) 
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